Topological insight into the non-Arrhenius mode hopping of semiconductor ring lasers 
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We investigate both theoretically and experimentally the stochastic switching between two 
counter-propagating lasing modes of a semiconductor ring laser. Experimentally, the residence time 
distribution cannot be described by a simple one parameter Arrhenius exponential law and reveals 
the presence of two different mode-hop scenarios with distinct time scales. In order to elucidate 
the origin of these two time scales, we propose a topological approach based on a two-dimensional 
dynamical system. 

PACS numbers: 42.65.-k,42.55.Px,42.60.Mi 



Fluctuations in active optical systems such as lasers 
is one of today's technological challenges as well as a 
fundamental problem of modern physics as they are the 
result of the quantum nature of the interaction between 
light and matter [TJ. Fluctuations are e.g. responsible 
for longitudinal mode switching in semiconductor lasers 
[2], polarization mode-hopping in Vertical Cavity Sur- 
face Emitting Lasers (VCSELs) [3H5], and they play a 
fundamental role in stochastic and coherence resonances 
of optical systems [SHE]- 

Semiconductor ring lasers (SRLs) are a particular class 
of lasers whose operation is strongly affected by stochas- 
tic fluctuations. The circular geometry of the active cav- 
ity allows a SRL to operate in two possible directions, 
namely clockwise mode (CW) and counter-clockwise 
mode (CCW). From the application point of view, SRLs 
are ideal candidates for all-optical information-storage. 
[9lfTT]. From a theoretical point of view, SRLs represent 
the optical prototype of nonlinear Z 2 -symmetric systems 
[12) , which appear in many fields of physics. 

Fluctuations induce spontaneous abrupt changes in the 
SRL's directional operation from CW to CCW and vice 
versa, and therefore represent a major limitation to their 
successful applications for instance as optical memories. 
An in-depth understanding of the mode-hopping in SRLs 
would shed light on the stochastic properties of the large 
class of ^-symmetric systems. In spite of its importance, 
the problem of spontaneous directional switches in SRLs 
remains unaddressed, partly due to the high dimension- 
ality of the models that have been proposed for SRLs 

mm- 

In this paper, we address the problem of such fluc- 
tuations both theoretically and experimentally. We ex- 
perimentally investigate the properties of the residence 
time distribution (RTD) that quantifies the mode hop- 
ping. Our theoretical analysis is based on an asymptotic 
reduction of a full rate-equation model to a ^-symmetric 



planar system [14] . 

We consider here an InP-based multiquantum-well 
SRL with a racetrack geometry and a free-spectral- 
range of 53.6 GHz. The device operates in a single- 
transverse, single-longitudinal mode regime at wave- 
length A = 1.56y«m. However, it will be clear from the 
rest of the discussion that our analysis is general and ap- 
plies to any kind of ring geometry. A waveguide has been 
integrated on the same chip in order to couple power out 
from the ring. This bus waveguide can be independently 
biased in order to reduce absorption losses. The waveg- 
uide crosses the facets of the chip under a 10° angle in 
order to minimize back-reflections. The chip is mounted 
on a copper mount and thermally controlled by a Peltier 
element which is stabilized with an accuracy of 0.01°C. 
The power emitted from the chip is coupled to a multi- 
mode fibre and detected with a 2.4GHz photodiode con- 
nected to an oscilloscope. We forward-bias the waveguide 
at a fixed current of 8.22 ± 0.01 mA in order to achieve 
transparency, and we bias the device with increasing DC 
current until it reaches the threshold at « 31.5 mA. The 
stochastic mode-hopping starts at approximately 39 mA. 
In the mode-hopping region, the time series of the power 
emitted in the CCW mode have been digitally recorded 
with the oscilloscope for different values of the bias cur- 
rent. Due to the anti-correlation of the two counterprop- 
agating modes [5] , a drop in power in the CW mode cor- 
responds to an increase in power emitted in the CCW 
mode. The CCW mode's RTD has been obtained from 
the time series and is plotted in Fig. [TJb). It is evident 
from Fig. |T|(b) that two well-separated time scales are 
present in the RTD. The white lines depict the best lin- 
ear fit in this logarithmic plot of the regions where the 
residence times are exponentially distributed with one of 
both characteristic time scales [cx exp (— i/(T 12 ))]. The 
fast time scale corresponds to an average residence time 
(Ti) of 10ns, whereas the slow time scale events have an 
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FIG. 1. (a) Schematic sketch of the set-up. (b) Mea- 
sured RTD in the CCW mode for a laser current of J,-i„ 9 = 
39.92mA; the waveguide is biased at 8.22mA and the temper- 
ature is 21.55°C. (c) RTD obtained by simulating Eqs. (|TJ- ([2]) 
with the following parameters: /j, = 1.59, a = 3.5, s — 0.005, 
c = 0.01, k = 0.44ns" 1 , fa = 1.5, D = 6.5 ■ 10~ 5 ns~\ Best 
fitting lines are shown in white. 

average residence time (T 2 ) of approximately Ifis. An 
increase in the bias current affects the time scales in dif- 
ferent ways: (Ti) remains almost constant, whereas (T 2 ) 
increases with bias current. The appearence of the short 
time scale is a characteristic feature of SRLs: similar ex- 
periments performed in VCSELs [4] [5] or in Dye ring 
lasers [16] lead to RTDs that are well fitted by a one- 
parameter Arrhenius distribution. A spectral analysis of 
the time series reveals a relaxation-oscillation (RO) fre- 
quency between 1.0 and 2.0 GHz and therefore rules out 
ROs as the origin of the non-Arrhenius features of RTD. 

In order to clarify the origin of the two separate time 
scales, we have directly investigated the recorded time 
traces. Qualitatively different switches are revealed in 
these time traces as shown in Fig. [^a)-(c). Fig. ga) 
shows the system leaving the initially lasing CW mode 
and settling into the CCW mode. Short excursions of 
the duration of approximately 10ns to the opposite mode 
have also been observed as shown in Fig. [2jb) . Opposite 
to the case of Fig.[2ja), during these short excursions, the 
system does not settle in the other mode. The durations 
of such events match the fast time scale (Ti) in the RTD. 
Instead, the second time scale in the RTD is related to 
long residences in one mode such as depicted in Fig.[2ja). 
The long residences are distributed in an exponential way 
with a mean-residence-timc of 1.16/ns. We would like to 
stress here that the RTD is not a simple superposition 
of independent events such as those in Fig. [2ja) , (b) , but 
that more complex transitions are also present in the sys- 
tem, such as the one depicted in Fig.[2^c). In the latter, a 
number of consecutive, correlated short time excursions 
are performed before the system settles in the counter- 
propagating mode. 

In order to explain our experimental results, we start 



from a classical rate-equation model, which accounts for 
both saturation and backscattering effects [9 3 . The evo- 
lution of the slowly varying envelopes of both complex 
counter-propagating fields (Ei and E 2 ) and the carrier 
density (N) is described by the following set of equations: 

£1,2 = k(1 + ia) [N (1 - s|£i, 2 | 2 - c\E 2 .i\ 2 ) - l] E h2 
-ke^«E 2A , (1) 
N = 1 [»-N-N(l-s\E 1 \ 2 -c\E 2 \ 2 ) \E^ 2 

-^(l-si^i 2 -^] 2 )]^! 2 ], (2) 

where the electric fields and the carrier inversion have 
been non-dimcnsionalized and the dot represents differ- 
entiation with respect to time. In Eqs. ([lJ-Q, k is the 
field decay rate, and 7 is the decay rate of the carrier 
population, a is the linewidth enhancement factor, fi 
is the renormalized injection current (ft w at trans- 
parency, /j « 1 at lasing threshold). Self- and cross- 
saturation effects are the result of spatial and spectral 
hole burning in the system and are described by s and c. 
The backscattering (for instance due to defects or back- 
reflections |151 ) is characterised by a certain amplitude k 
and phase (f>k . Although perfectly adequate to model the 
SRL operation, nevertheless such five-dimensional model 
needs to be simplified in order to allow for a clear topo- 
logical analysis. A reduced two-dimensional model has 
been proposed [14], which describes the asymptotic be- 
haviour of the ring laser on time scales that are slower 
than the RO: 

9 = — 2 sin (fib sin ip + 2 cos 4>k cos ?p sin 9 

+ J sin 9 cos 9, (3) 
cos 9ip = aJ sin 9 cos 9 + 2cos(f>k sin ip 

+ 2 sin 4>k cos ip sin 9. (4) 

where 9 = 2 arctan ^paWPiF - tt/2 € [-tt/2, tt/2] 
represents the relative modal intensity and ip € [0, 27r] 
is the phase difference between the counter-propagating 
modes. The pump current has been rescaled as J = 
«(c - s)(n - I) Ik. The use of Eqs. Q-Q is justified by 
the fact that all the observed time scales [Fig. [ljb)] are 
slower than the relaxation oscillations. 

When the SRL operates in a unidirectional regime, 
the phase portrait of the system is exemplified in Fig- 
ure |2^g)-(i). Four stationary solutions exist for Eqs. ([3])- 
Q: an unstable in-phase bidirectional state in (0, 0) (not 
shown); two symmetric stable states CW and CCW at 
ip 7r, both corresponding to unidirectional operation; 
and a saddle point S in (0, it) which is the unstable out- 
of-phase bidirectional solution. The stable manifold of 
S separates the basins of attractions of CW and CCW, 
whereas the unstable manifold connects S with CW and 
CCW. At a critical value Jhom of the current, a ho- 
moclinic bifurcation takes place in the system and the 
stable and unstable manifold of S coincide [14]. When 
J > Jhom the stable manifold spirals around S and the 
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FIG. 2. Comparison between measured time series (a)-(c), simulated time series (d)-(f) and phase space trajectories (g)-(i) 
for different kind of transitions. In the experiment the ring was biased at J r i ng = 39.86mA; the waveguide was biased at 
8.22mA and the temperature was 21.55°C. In the numerical simulation the following parameters were used: fi = 1.59, a — 3.5, 
s = 0.005, c = 0.01 k = 0.44ns -1 , cj> k = 1.5, D = 6.5 ■ 10" 5 ns _1 . The notation of (g)-(i) is as follow: CW, CCW are the stable 
states, S is the saddle; solid line: stable manifold of S, dashed line: unstable manifold of S, gray: projection of the time series 
(d)-(f). 



basins of attraction of CW and CCW fold around each- 
other. Increasing the parameter J continuously unfolds 
the stable manifold. 

When noise is present in the system, a rare, large fluc- 
tuation may drag the system away from an initial sta- 
tionary state to the basin of attraction of the opposite 
mode. It is known from the theory of stochastic transi- 
tions in nonlinear systems that in the limit of small noise 
strength a transition takes place in a ballistic way along 
a most probable escape path (MPEP) which can be cal- 
culated by solving an auxiliary Hamiltonian system 17 \. 
As our arguments rely only on the topological features of 
the system, an exact calculation of the MPEP is not re- 
quired. In a system such as Eqs. Q-Q, it is known that 
the MPEP connects a stationary state with the saddle 
S [18]. More specifically, the MPEP approaches S along 
its stable manifold. Once the saddle is reached, the tran- 
sition to the opposite equilibrium state is completed by 
following the unstable manifold of the saddle. However, 
when the system's parameters are close to the homoclinic 
bifurcation, the folds of the stable manifold cluster very 
tightly around each other, and the stable and unstable 
manifolds of the saddle become very close. Therefore, 
even at very low noise intensities, the system can diffuse 
between the folds of the stable manifold or between the 
stable and the unstable manifold. Such topological con- 
siderations can explain the separated time scales (Ti^) as 



well as the features of the transitions shown in Fig. [2ja)- 
(c) . A noise induced activation is responsible for the slow 
time scale (Ti), whereas a noise-sustained rotation along 
the folds of the stable manifold of S produces the short 
excursion shown in Fig.|2^b). 

We stress that SRLs are out-of-equilibrium systems 
[TT?] ; although the phase portrait of the system is sim- 
ilar to the one obtained in the case of an underdamped 
Duffing oscillator, the topology of the manifold is related 
to the vicinity of the homoclinic bifurcation, but cannot 
be considered as a sign of weak dissipation. 

In order to confirm our topological arguments, we per- 
form numerical simulations of Eqs. ([T|)-(|2| using relevant 
parameters [3] . The spontaneous emission noise has been 
introduced in Eq. (JlJ phenomenologically as complex un- 
corrected zero-mean stochastic terms described by the 
correlation terms: (£j (t + r) £* (t)) = 2DSijS (t) where 
i,j = 1,2 and D is the noise intensity. Carrier noise has 
been disregarded as its relevance in directional mode- 
hopping was proven negligible |16j . 

Three examples of numerically obtained stochastic 
transitions are reported in Fig. [2^d)-(f) and compared 
with the experimental ones [? ]. The same trajectories 
projected in the 9 — ip plane are shown in Fig. ptg)-(i) 
together with the invariant manifolds. The trajectory in 
Fig. [2^d) shows a simple mode-hop as in Fig. [2ja). In 
phase space [Fig. [2^g)] such trajectory corresponds to a 
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stochastic path which approaches the saddle and relaxes 
to the opposite state along the unstable manifold of S. 
A short excursion to the opposite direction is shown in 
Fig. [2|e), corresponding to the experimental time trace 
shown in Fig. [2^b). In phase space [Fig. [2jh)] the tra- 
jectory rotates around the saddle, but remains inside the 
basin of attraction of the initial mode. In Fig. |2jf), we 
show an example of a more complicated transition sim- 
ilar to the experimental trajectory shown in Fig. [2^c). 
In phase space [Fig. |2ji)] the system rotates three times 
around S before settling in the opposite state. The sim- 
ulations shown in Fig. |2^d)-(f) are in good qualitative 
agreement with the experiments in Fig. [2ja)-(c). A di- 
rect simulation of the reduced model Eqs. leads 
to similar results. 

We remark that the described effect takes place in re- 
gions of the parameter space (J,a,<f>k) which are close 
to the homoclinic bifurcation. In the unphysical case of 
k = (no backscattering) the manifolds are completely 
unfolded for every value of 4n- and J and the effect dis- 
appears. 

Finally, we have extracted the RTD from the simulated 
time series [Fig. [ljc)]. The simulations show the same 
non-Arrhenius structure as in the experimental RTD, and 
quantitatively reproduce the time scales (Ti^)- 

In conclusion, we have shown that the RTD for the 
mode-hopping in SRLs is strongly non-Arrhenius and 
cannot be described by a single transition rate. Two sep- 
arate time scales are observed in the RTD. Our explana- 
tion of this behaviour is based on topological properties 
of the phase space of the asymptotically reduced model 
[Eqs. ([Sjl-Q] in the vicinity of a homoclinic bifurcation. 
We have found that during a noise induced transition, the 
system explores the region of the phase space where the 
stable manifold of the saddle S rotates around the sad- 
dle. Even at low noise intensity, the system can diffuse 
between different folds of the manifold leading to different 
kind of transitions such as those experimentally observed 
in Fig. [2|a)-(c). Numerical simulations agree well with 
the experiments. The simulated RTD is in good quanti- 
tative agreement with the experimental one. We finally 
remark that the topology of the manifolds that we have 
observed for the planar system [Eqs. Q-Q] has been 
observed in other optical system such as in the Ikeda 
map applied to a nonlinear ring resonator |20j , and it is 
a general feature of two-dimensional dynamical systems 
which are Z 2 -invariant [IS]. Therefore, SRLs are excel- 
lent optical prototypes for the experimental exploration 
of stochastic systems with ^-symmetry. 
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